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' Abstract. Light storage in an atomic Bose-Einstein condensate is one of the most 

practical usage of these coherent atom-optical systems. In order to make them even 
more practical, it is necessary to enhance our ability to inject multiple pulses into 
' the condensate. In this paper, we report that dispersion of pulses injected into 

the condensate can be compensated by optical nonlinearity. In addition, we will 
present a brief review of our earlier results in which enhancement of light storage 
^ ' capacity is accomplished by utilizing multi-mode light propagation or choosing an 

optimal set of experimental parameters. 

p ; 

^' 

1 Introduction 

(N 

[ Soon after the generation of Bose Einstein condensate in an ultracold gas of trapped Alkali 

CS| i atoms [T], a promising utilization of it for the dramatic slowing down of a light pulse was 

demonstrated [2]. This feat is accomplished with the help of a quantum coherent effect called 
electromagnetically induced transparency (EIT) |3|4| . Ultraslow light pulses can be used for 
storage of coherent optical information [5J . Mutual conversion of classical coherent information 
I (phase and amplitude) of the light pulse and the quantum information (quantum state) of the 

atomic system can allow for quantum information processing via ultraslow light [6]. A quite 
recent experiment provides strong hope towards this direction [7j. In the experiment, coherent 
optical information is first encoded in one condensate. Quantum information in the condensate 
is then carried to another condensate by a matter wave. Finally, coherent optical information 
is revived in an optical pulse, generated out of the new condensate on demand. 
■ Due to such impressive developments in experimental ability to control of light and matter 

' waves in light storage experiments, more practical quantum information processing applications 

can be expected to occur in near future. On the other hand, it is necessary first of all, to increase 
our capability to control the amount of information stored in the condensate. For that aim, we 
should be able to inject more than one pulse into the condensate during the storage time. This 
is a basic requirement to realize practical logic gates in atomic condensates using slow light set 
ups |5|6j . 

To investigate how to make more than one pulse simultaneously present in an atomic con- 
densate efficiently, we have performed a series of studies. This paper reviews some of our earlier 
results and in addition it reports our new results where we have found optimum conditions for 
dispersion compensation via nonlinearity. In our earlier works, we show that optical dynamic 
memory capacity of the condensate can be optimized by choosing a certain set of experimental 
control parameters, such as coupling laser Rabi frequency and temporal width of the probe pulse 
for a given atomic condensate |£j . It has been shown that axial density profile of the condensate 
helps to preserve the probe pulse shape against group velocity dispersion [8] . Further enhance- 
ment of the memory capacity can be possible by taking into account radial confinement of the 
probe pulse. We have demonstrated that a particular radial density profile of the condensate 
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reduces the effect of modal dispersion and contributes to the pulse shape preservation [9]. We 
deduce that the optical control parameters are beneficial for optimizing the memory capacity 
and at the same time, properties of the atomic cloud can also be exploited to enhance the maxi- 
mum capacity available. When we take into account radial confinement, ultraslow wave-guiding 
regimes have been investigated and characterized. Propagation constants and the conditions on 
the number of ultraslow optical modes have been determined. Detailed numerical examinations 
revealed the ultraslow mode profiles that can be supported by atomic condensates. In this pa- 
per, we will go beyond these linear optical results and examine the effects of nonlinear optical 
response of the condensate. We show that ultraslow solitons can bring further enhancement in 
the light storage capacity. 

Following this introduction, in Sec. 1, the paper first reviews optical information storage 
enhancement schemes based upon linear optics considerations. In Sec. 2.1., group velocity dis- 
persion and a scheme to beat it are discussed. In Sec. 2.2, taking into account transverse 
directions to the propagation direction, modal dispersion is studied. Information storage en- 
hancement with the aid of ultraslow waveguiding modes are discussed. Our new contributions 
in which ultraslow light propagation is studied in the nonlinear regime is described in Sec. 3. 
Increase in the capacity of optical information storage due to the ultraslow optical solitons is 
discussed in that section. Finally, we conclude in Sec. 4. 



2 Control of optical information storage in linear optical regime 

To slow down group velocity of an optical pulse just to few meter per second, steep dispersion of 
electromagnetically induced transparency is used [2] . Optical pulse to be slowed down consists 
of a weak probe light. It is shined upon atomic BEG, together with a stronger control field. This 
configuration allows for cancellation of absorption at resonance and creates transparency over a 
small frequency window. Measuring the axial length of the condensed cloud by an imaging laser 
and comparing the time delay in the probe beam propagation with a reference beam not entered 
into the condensate determines operationally the average group speed of the probe pulse. Such 
a typical slow light experiment is usually modelled by a one dimensional linear optics model. 
This model predicts the measured group speeds well, especially when the interactions within 
the cloud is taken into account ^lOj . 

Considering a gas of effectively three-level atoms, linear optical response of an FIT medium 
is characterized by a susceptibility given by 

_ pM i{iA + r^/2) 

^ eoh [{r2/2 + iA){r3/2 + iA} + n^j4y ^' 

The effective level configuration of the atoms consists of two dipole allowed transitions such 
that the lower levels |1) and |2) are coupled to the upper level |3) by the probe and the control 
fields, respectively. Here, p is atomic (number) density. We assume that the control field of 
Rabi frequency f2c is at resonance, while A is the detuning of the probe field of frequency uj 
from the resonance frequency luq- M3i is the dipole matrix element for the probe transition. 
Taking resonance wavelength A for the probe transition we have the usual relation p.^i — 
3eohX^j/8Tr^ with 7 is the spontaneous decay rate for the probe transition. Dephasing rates of 
the corresponding transitions are denoted by / 2 and la . Within a small frequency window about 
the probe resonance, susceptibility becomes a negligibly small pure imaginary number, which 
leads to the transparency effect. Besides, within this window, susceptibility shows steep variation 
with A so that group velocity reduces dramatically for a resonant probe pulse that can propagate 
without suffering any significant absorption. Physically, this is possible by the cancellation of 
absorption due to the destructive quantum interference of the transition probabilities to the 
upper level. Quantum coherence is established by the control field which is stronger than the 
probe field in the usual FIT condition. 

It may be argued that BFC is not essential for the ultraslow light generation. Degenerate 
Fermi gas also allows for significant reduction in group velocity, though not as efficient as BEG 
[TT] . BEG gives much slower speeds over more compact spatial dimensions. Its unique advantage 
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however is the abihty to combine matter wave dynamics with the uhraslow light pulses [7] . This 
seems to be more promising for practical quantum information processing. 

Susceptibility formula is derived by the steady state solution of the Liouville equation for 
the density matrix, assuming most of the atoms remain in level In the subsequent section, 
we will go beyond linear regime and consider nonlinear response. In the linear regime, we will 
first discuss dispersive propagation of ultraslow short pulses in one dimension. After that, we 
extend this by taking into account transverse directions as well. 

Dispersive effects are weak in the usual ultraslow light experiments as these pulses are just 
a few microseconds long. First attempts to model experimental system uses one dimensional 
non-interacting BEC as the propagating medium with linear optical response [12] . Further- 
more it is assumed that the medium is non-dispersive. This model is improved significantly 
by considering atom-atom interactions. It has been shown that increasing interactions lead to 
faster group speeds |T^ . Higher dimensional propagation effects within paraxial approximation 
are also considered and diffraction effects at the thermal cloud boundary of the condensate 
are numerically simulated |10| . These results demonstrate validity of the non-dispersive, one- 
dimensional model in describing experimental results. On the other hand, when it comes to 
improve storage capacity, this model becomes insufficient. 

Natural attempt to enhance our ability to inject multiple pulses into the BEC during the 
storage time is to use shorter pulses in time. By this way, if we say every pulse carries a bit of 
coherent optical information (which is the classical information associated with the amplitude 
and phase of the pulse), then we could achive higher bit storage capacity. With the usual 
microsecond pulses however, as the storage time is also about a few microseconds, one can get 
only one pulse present in the condensate [2] . Some schemes are proposed to increase number of 
pulses. One can utilize for instance different polarization of probe pulses or enhancement of 
EIT window [T3] . Both these proposals would benefit from more direct approach of being able to 
inject shorter pulses into the condensate |15|16|17] . The fundamental problem in using shorter 
pulses is that they may be broadened due to the group velocity dispersion within the storage 
time. It is thus necessary to study effect of group velocity dispersion in the axial propagation 
of ultraslow short pulses in order to estimate actual enhancement of the bit storage capacity. 



2.1 Group velocity dispersion 

We will now investigate the propagation of short laser pulses under slowly varying phase and 
envelope approximations without taking into account nonlinear optical response. A formulation 
of light propagation through a dispersive EIT medium is given in Ref . |18| . We have applied this 
theory to ultraslow short light pulses through BEC in Ref. |8j . We now review and summarize 
our major results in Ref. [8J. 

Total polarization is related to the electric field through electric susceptibility via P = P*^') 
where the linear polarization is P*-^^ = eox('-^)^^^E. Dispersive medium is characterized by a 
frequency dependent susceptibility. We write P{lj—ljo, E) ~ eoxl^^— ^o)^^^-E'('^— wq). Expanding 
the linear susceptibility x^^^ of f^ic dressed atom to second order about a central frequency 

polarization becomes 

P{t) = .ox(-o)( ^E{t) - ^eo^U^ - J^I^U^- (3) 



In slowly varying approximation, the wave-equation is given by. 
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Using equation ^ in Eq. HJ we obtain the wave equation for the short optical pulse such that 

dE ^ 1 dE d^E 

where a determines the pulse attenuation; Vg is the group velocity, 62 is the second-order dis- 
persion coefficient, or the group velocity dispersion. Calculating the third and the second order 
dispersion coefficients using typical experimental parameters [2] but for shorter pulses, we have 
found that third order dispersion is seven orders of magnitude smaller and negligible. Second 
order dispersion becomes significant for pulses shorter than microseconds. In our numerical 
calculations, we consider parameters for the experiment in Ref. 2 except for the coupling field 
Rabi frequency which is chosen to be Qc — 87. The transmission window, Alu is calculated to 
be Au! = 8.49 x 10^ Hz so that we shall consider pulses of width up to r = 10 ns for efficient 
transmission. The coefficients can be calculated directly from the susceptibility via 

ZTT 

a = -yxC^^o), (6) 



1 TT 9x 
c A duj 



^2 - — ITTjI'^o- 



These expressions lead to complex parameters in general, as x is a complex valued function. 
At EIT resonance, they lead to physically meaningful well-defined real valued absorption coef- 
ficient, group velocity and dispersion coefficient. When Rabi frequency for the coupling field is 
sufficiently large such that f2c ^ -^2,3, the coefficients reduce to 

^ 2tJ3i|^3l|2p' 

^ . 87rr3|M3i|V 

It should be noted that the given expression for for Vg is valid only when 



(10) 

(11) 
(12) 



r2,3 < < ^JA^Tc\^l:il\^^p|eon = ^cX^lP- (13) 

From the explicit expressions, we see that Vg increases more slowly with the Rabi frequency of 
the control field Qc than the decrease of the dispersion coefficient 62 . This higher susceptibility 
of 62 to Qc can be exploited to optimize the bit storage capacity C, defined by C = L/2vgT 
with T is the largest width of the probe pulse within the condensate of size L. 

Analytical solutions for C are derived for a Gaussian pulse assuming a uniform effective 
dispersive medium representing the actual spatially inhomogeneous BEC [8, . Analytical solution 
of the wave equation allows for determination of the pulse width. L is calculated as the root- 
mean-square of the axial distance for a given condensate of number density p, under a semi-ideal 
model of the weakly interacting condensate . This leads to 

C= ^ (14) 

2to V4^'^^c/9A472p2 + AL^r'i/n-^T^Qt 

Initial pulse width is denoted by tq . C shows a maximum at a critical Rabi frequency determined 

by 

^co = ^2^2 • (15) 
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The corresponding storage time of the pulse t^Q is given by 




From the number of pulses that can be present in this storage time we find the maximum 
capacity as, 

^vaax - y 32^^ , 

which is independent of the pulse width at the critical Rabi frequency. At the critical Rabi 
frequency, about two orders of magnitude improvement is achieved regarding the bit storage 
capacity. 

More rigorous numerical calculations are verified these analytical results. Numerical calcu- 
lations take into account density variations and the thermal component of the gaseous system. 
Temperature and the interactions affect the pulse character through the dependence of pulse pa- 
rameters on spatial size and density product. Numerical studies are required to handle spatially 
inhomogeneous BEC. Propagation parameters in the wave equation gain non-local character 
under this model but can be treated numerically. We have found that as the pulse gets broader 
and broader as it propagates towards to the central region of the atomic cloud, where the dis- 
persive effects are most strong, it becomes better protected against dispersion. As it is already 
broadened, it gets less broadening effect. Axial density profile of the condensate protects the 
pulse against group velocity dispersion. Thermal variation of the bit storage capacity is investi- 
gated numerically and just below critical temperature, an optimum temperature for the highest 
capacity is determined. Interactions in the cloud effects the broadening of the pulse and makes 
it more sharp when the interactions are stronger. Fig. [T] illustrates these results. 




Fig. 1. Temperature dependence of the pulse width calculated numerically for an inhomogeneous 
condensate profile. Initial width of the pulse is 1 ns. Effect of the interactions can be seen by comparing 
the solid curve for which the s-wave scattering length characterizing the effective atom-atom interactions 
tts = 2.75 nm with the dashed one for which = 7 nm. The parameters used in the calculation are 
taken from the experiment in Ref. i2j. 

From the figure we see that the short pulses suffer significant broadening despite the short 
size of the condensate. One cannot increase the capacity C thousand times by using nanosecond 
pulses instead of microsecond ones. Using a control field at the suggested critical value as above, 
an optimum increase of the number of pulses that can be injected multiply into the condensate 
can be achieved. Analytical result of C evaluated gives about 80 pulses in the condensate while 
numerical result suggests, as explained above, due to less broadening within the inhomogeneous 
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condensate, 800 pulses could in fact be stored. In comparison to present slow light experiments, 
that means a two orders of magnitude enhancement in the bit storage capacity. Below we 
shall discuss additional improvements to C, taking into account transverse directions of the 
condensate as well as nonlinearity in the probe pulse propagation. 

2.2 Modal dispersion 

In a subsequent analysis we take into account transverse dimensions to the axial propagation 
direction and examine the ultraslow optical modes. Waveguiding problems in cold gaseous 
systems was subject to much attention recently |20I21I22I23| . Now we summarize our major 
results reported in Ref. [9]. We observe that at extremely low temperatures well below the 
critical temperature of condensation, density profile of the condensate can be approximately 
described by Thomas-Fermi approximation. Normally, at perfect EIT resonance, no significant 
index difference between the condensate core and the thermal background gas can be achieved. 
In this case, only single mode is supported in the condensate. For a slightly off-resonant probe 
however, sufficiently high index contrast can be generated between the condensed core and 
the enveloping thermal gas. Under these conditions, radial density of the condensate can be 
translated to a graded index profile characterized by a parabolic (quadratic) spatial variation 
in the form 



with A = 1 — 1 /nf. Here R is the thermal radius along the transverse (radial, r) direction and ni 
is the index along the axial (z) direction. Such a system resembles a graded index fiber, though 
one of micron sized. Light propagation is found to be in the weakly guided regime where the 
optical modes are described in terms of linearly polarized (LP) modes as the index difference 
is small and the axial component of the electric field is negligible. In this case the radial wave 
equation is described by the Helmholtz radial equation 



Here ko — uj/c,l — 0,1,2,..., /?, and ip are defined for the transverse field of the LP modes 
Et = tpir) exp + ujt - f3z))]. 

The number of such transverse (LP) modes as well as their group velocity is evaluated. 
We find that modal dispersion is stronger at low temperatures and large number of modes 
can be supported by the condensate. As the temperature gets higher and higher, number of 
modes that can be supported decrease. This is due to the shrinkage of the condensate and 
the vanishing of the index contrast. Analytical results are provided by a conventional WKB 
treatment. Conditions for multiple, single and two mode propagations are provided. It is shown 
that all LP modes can propagate at ultraslow light speeds. Lower order modes propagate 
slower while higher order modes goes faster as they are away from the central, more dense, 
condensed region. This particular density profile of the BEC then allows for less pulse spread 
due to differences in the group speeds of different modes. As in the case of axial propagation 
where axial density of the condensate protects the pulse against group velocity dispersion, 
now the radial profile of the condensate also helps for the pulse shape preservation against 
modal (waveguide) dispersion. More detailed numerical studies are performed. We consider the 
condensate in terms of tiny shells and treat the index as constant within each region between 
the shells. At constant index, analytical solutions in terms of Bessel function are written down 
and matched at the boundaries. By this way, actual mode profiles are numerically determined. 
We illustrate two modes in Fig. [51 

Short axial size of BEC and its unique density profile can be exploited to keep modal 
dispersion weak. In such short distance applications, this may be utilized to capture optimum 
amount of power from a light source and enhance the optical coupling/sensing capability of 




(18) 



(19) 



BEC. 
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r/R 



Fig. 2. Two ultraslow transverse modes confined in a BEC. Solid curve is for LPqo while the dotted 
one is for LPio. Mode profiles are calculated with the parameters that are taken from the experiment 
in Ref. [2], except for A = — O.I7 and Sic — 2.67. 



LP modes can be addressed by specifically constructed transverse profiles of tlie probe 
pulses. This can be exploited to control of capacity in the transverse direction for storage of 
coherent optical information in different mode patterns. These results can be further combined 
with the longitudinal control of optical information storage capacity where a particular EIT 
scheme, in which the control field is perpendicular to the probe field, is employed [24j . 



3 Control of optical information storage in nonlinear optical regime 

We will now take into account nonlinear optical response of the condensate in the short 
pulse propagation. Total polarization, now including the contribution of nonlinear polariza- 
tion P("') = eo(x(a;)(2' : EE + xi^)^^^ ■ EEE), can be expressed as 

P{lu - LOO, E) = eox(t^ - loqY^'^E{lo - ua) + eox^^^ : EE + e^x^^^ : EEE, (20) 

where an expression for the third order nonlinear susceptibility is given in Ref. |25| . This 
introduces nonlinear terms to the susceptibility expansion so that polarization can be written 
as 

P{t) = .ox(-o)^ ^E{t) - ^eo^U^ - -^^U — 

+ eox'^' : EE + eox'''^' : EEE. (21) 

X^"^^ is zero due to the inversion symmetry of BEC. Using Eq. (PT|) in Eq. (H]), we get a nonlinear 
wave equation for the short optical pulse such that 

dE 1 dE d^E , ,„ 

^+-E+--+^b.— +^,\E\^E = 0, (22) 
where 77 is the nonlinear optical Kerr coefficient that can be calculated from the susceptibility 



via 



When i7c > 3, we find 



ry = --x(^)(c.o). (23) 



_ 4^p|/i3l|'/^2 (1^2 + ^3) 
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Let us now write the wave equation ([22]) in a more conventional form |26j , 

Here /3i — P2 — 262 and a —2a. By making transformation T = t — z/vg, A = ^t^cE 

and 7 = Tj/tQC [26j we obtain 

^a^ + ^T^-2^^^-^l^l^ = «- (2^) 

where A is the slowly varying amplitude of the pulse envelope. 

Spatial variation of the nonlinear optical Kerr coefficient in the axial direction is shown in 
Fig. ([3]). Fig. (|4]) displays the axial spatial profile of the second-order dispersion coefficient. 
These figures demonstrate a potential for compensation of the dispersion by nonlinearity is 




-200-100 100 200 
z(jum) 

Fig. 3. Axial spatial profile of the nonlinear coefficient 7 for the parameters same with those of Fig[T] 
7 is normalized by its peak value 3.64 x 10~'' m/W. 



possible. We need to determine required laser power for that aim. 




-200-100 100 200 
z(yum) 

Fig. 4. Axial spatial profile of the second-order dispersion coefficient P2 for the parameters same with 
those of Fig[T] /?2 is normalized by its peak value 1.07 x 10"^^ s^/m. 



In order to solve Eq. (|26p one can introduce normalized amplitude U by using the definition 

m 

A{z,t) y/Poexp{-a' z/2)U{z,t). (27) 

where Pq is the peak intensity of the incident pulse and r — T/Tq with initial pulse width 
Tq. The dispersion length Ljy and the nonlinear length L^l provide length scales over which 
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dispersive or nonlinear effects become important for the pulse evolution in the Bose-Einstein 
condensate with an effective length L. The dispersion length Ljj and the nonlinear length L^l 
are given by [26] 

L. ^ ||, (28) 
For the nonlinear pulse propagation, a dimensionless parameter N is introduced as 

^'-T^-'-^- (30) 

i^NL \P2\ 

For values TV « 1 both Self-phase modulation(SPM) and Group-velocity dispersion(GVD) play 
an equal role [26] . Under this condition, we determine the peak power to compensate the group- 
velocity dispersion for a short optical pulse propagating in the BEG using 



^0 = ^. (31) 



Using Eqs. pT]l and ([24]) . we find the relation between the probe pulse characteristics and the 
material parameters such that 

(32, 

Using the parameters given in the slow light experiment Ref. such that 7/27r = 10.01 MHz, 
A — 589.76 nm, J3 = 7/2, r2/27r = 1 kHz, we calculate that Pq — 389 mW/cm^ for a 1 /iS pulse. 
For a shorter pulse of 10 ns, required intensity increases to 3.9 kW/ /cm^. These intensities are 
beyond the saturation intensity of the corresponding transition for the sodium atom, which is 
in the order of ~ 10 mW/cm^. However, it is common to consider high intensities, in particular 
for the control field in EIT. For example, Ref. [25] considers a control field of intensity 55 
mW/cm^. In the EIT process, absorption of the probe field is cancelled due to the quantum 
interference of the transition amplitudes. The major limitation on the probe intensity is the 
EIT condition itself, which requires much higher control field Rabi frequencies than that of 
the probe field. If we use the same parameters with Ref. [25], which takes slightly different 
decoherence and decay rates than the ones used in the experiment in Ref. [5], we would get 
Po = 1.6 mW/cm^ for a 1 /zs pulse and Pq = 166 mW/cm^ for a 0.1 /is pulse. Sensitivity of the 
required peak intensity on the natural linewidths and the decoherence rate of the atomic system 
suggest that one can engineer linewidths using different atoms, level structures [30 and putting 
the atoms in an optical cavity [31j to optimize demanded peak intensity. Furthermore, instead 
of a continuous wave (cw) shining upon the condensate, we consider storage of pulse trains 
inside the condensate, and hence the average power would be smaller than the one associated 
with the peak intensity. The results reported here can be extended to the ultraslow pulses in 
hot gases as well, by taking into account Doppler and collision broadenings. In such systems, 
much higher pulse intensities in the orders of 1 — 1000 W/cm^ are used in EIT experiments 
[32]. 

We would like to emphasize that other applications based upon nonlinear effects might have 
different and less restrictive requirements on the peak intensity or power of the probe pulse. As 
a low loss medium with enhanced nonlinearity [33], with small spot size due to small conden- 
sate cross section, self-focusing effects and EIT on the probe beam focusing, many nonlinear 
effects of various interest could be realized with low powers in the condensates in the EIT 
configuration. Our focus and the purpose here is to examine the standard method of dispersion 
compensation via nonlinearity so that number of pulses that can be simultaneously present 
in the condensate can be optimized more directly. In the earlier sections, we have discussed 
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indirect methods based upon adjusting control field intensity as well as taking advantage of 
higher dimensional pulse propagation and the transverse optical modes inside the condensate 
to enhance the optical storage capacity. The present section demonstrates the possibility of 
exploiting nonlinear compensation. With less demanding probe pulse power requirements, rich 
physical situations can be created such as quantum correlations in the weak probe induced by 
the nonlinearity |25| . From classical point of view the relatively simpler theory of nonlinear 
pulse propagation is employed here. Beyond dispersion compensation, such a simple and well 
developed theory can also lead to rich applications and some novel effects, in particular in the 
context of dark and bright ultraslow soliton propagation [27128129] , coupled matter wave and 
optical field nonlinear dynamics, and so on. However their discussion is beyond the scope of the 
present paper, which is limited to optical storage capacity enhancement in BECs. 



4 Conclusion 

We have investigated the problem of injecting multiple ultraslow pulses through a BEC. We 
have examined three methods for that aim. This paper first reviews our two earlier results 
where nonlinear optical response is ignored, then also presents our new results where nonlin- 
earity is present. In our first study, we have considered short ultraslow pulse propagation in 
one dimension. Taking into account group velocity dispersion, an optimum set of experimen- 
tal parameters, in particular Rabi frequency of control field, that would give maximum pulse 
storage capacity are determined. Secondly, transverse direction is taken into consideration, and 
waveguiding of ultraslow pulses is studied. Existence of ultraslow transverse modes is revealed 
and their benefits in enhancement of pulse storage capacity is pointed out. Finally, in our new 
results, we have introduced nonlinear optical response and showed that it is possible to compen- 
sate dispersion of a short ultraslow pulse for a specifically chosen peak power of the pulse. This 
may be more efficient in optical storage than group velocity dispersion by the Rabi frequency of 
the control field. We hope these studies complement ongoing experimental efforts for practical 
storage of coherent optical information in atomic Bose-Einstein condensates. 
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